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Is the notion of a quantum computer resilient to thermal noise unphysical? We address this 
question from a constructive perspective and show that local quantum Hamiltonian models provide 
self-correcting quantum computers. To this end, we first give a sufficient condition on the connect- 
edness of excitations for a stabilizer code model to be a self-correcting quantum memory. We then 
study the two main examples of topological stabilizer codes in arbitrary dimensions and establish 
their self-correcting capabilities. Also, we address the transversality properties of topological color 
codes, showing that 6D color codes provide a self-correcting model that allows the transversal and 
local implementation of a universal set of operations in seven spatial dimensions. Finally, we give a 
procedure to initialize such quantum memories at finite temperature. 
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I. INTRODUCTION 



The issue of stability of self-correcting quantum mem- 
ories [H Q at finite temperature has recently attracted 
much attention. One reason for this is that they are more 
physical than the models of fault-tolerant quantum com- 
putation based on the quantum circuit model While 
the latter is based solely on the naked quantum formal- 
ism with no dynamics, the topological memory relies on a 
suitable Hamiltonian that protects quantum information 
as a classical information is protected in, e.g., a 2D ferro- 
magnet. The original topological memories are placed on 
a torus, and it is known that they offer stable qubits only 
starting from four spatial dimensions 0, 0j S B 0] ■ This 
may sound pesimistic, but there is a hope that the 3D Ba- 
con model [g| will work 0] . This model is not topological, 
but inherits the Hamiltonian approach from topological 
models. Yet, the Hamiltonian for Bacon model is not 
analytically diagonalisablc, and therefore its analysis is 
quite challenging. 

Both the Bacon and the toric models have some sim- 
ilar drawbacks. They do not support a set of universal 
gates, and no scheme is known for the initialization of 
these self-correcting memories in the Hamiltonian spirit. 
Namely, one would like to have something similar to an 
Ising model, where the state of definite magnetization is 
prepared by simply switching on a global magnetic field 
and then switching it off adiabatically. Thus, disregard- 
ing for a while the dimension of the model, an interesting 
question arises: does it exist a self-correcting quantum 
computer, allowing universal computation as well as a 
natural initialization scheme? 

In this paper we propose a first such model. To this 
end we use color codes, a class of topological codes with 
very special transversality properties [1, [TO] . Color codes 
are obtained from colexes, a class of lattices introduced 
in and their underlying mathematical structure is a 
simplicial homology that differs from the usual one be- 
cause it is 'colored'. In the color code Hamiltonian mod- 
els excitations take the form of branyons [ll|, punctual 



or extended objects subject to topological interactions. 
These models are local but require interactions between 
many spin-1 /2 s ystems, a difficulty that can be overcome. 
E.g., in Ref [l2] present an experimental proposal to im- 
plement the many-body interactions present both in toric 
and color codes. Alternatively, a local two-body Hamil- 
tonian model that effectively yields a two dimensional 
color code has been described recently [ist . 

The main idea of the paper is that in a suitable spatial 
dimension it is possible to find color codes that support 
universal quantum computing and, at the same time, pro- 
vide self-correction (which we prove giving a Peierls ar- 
gument in the style of We also present a scheme of 
initialization, analogous to that of a ferromagnetic mem- 
ory. We believe that this is the first proposal of such an 
initialization, as opposed to standard methods in fault 
tolerance (usually based on measuring syndrome on some 
initial product state). There is a remarkable difference 
between our initialization scheme and that of a ferromag- 
netic memory. Namely, we need to switch off at some 
point part of the protecting Hamiltonian, which is not 
needed in in the latter case. Interestingly, this is com- 
patible with a recent result of Alicki , who claims that 
if we were able to perform quantum computation with the 
protecting Hamiltonian always on, then we would violate 
the Second Law. 

The paper is organized as follows. Wc first recall 
the previous results regarding thermal stability of self- 
correcting quantum memories (Sec. |lT]). In Sec. IIIII we 
discuss the subsystem picture, which is natural for our 
analysis of thermal stability. In Sec. IIVI we describe 
general properties of stabilizer codes. Wc also introduce 
'dressed' observablcs - the ones that define the protected 
qubit. In Sec. |V]we present our modeling of the interac- 
tion with the environment. We recall an upper bound on 
the decay rate of the autocorrelation functions, and its 
relation with fidelity. Wc then formulate sufficient condi- 
tions for self-correction, related to the Peierls argument 
in Sec. IVII In addition, we prove that the given condi- 
tions also ensure that the quantum memory is resilient to 
a certain amount of noise applied suddenly. Next we turn 
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to the construction of the models, which must satisfy the 
conditions and at the same time allow a universal set 
of transversal operations (Sec. IVIip . We give a unified 
homological picture for toric codes and color codes, and 
prove in detail that some color codes possess both prop- 
erties at the same time. Finally, in Sec. I Villi we propose 
a scheme for initialization, measurement and gate appli- 
cation, arguing heuristically why it should work. The 
scheme relies on the ability to switch off the protecting 
Hamiltonian for a short time. 



II. PREVIOUS RESULTS 

The idea of a thermally stable self-correcting quantum 
memory goes back to The authors propose a 4D 
topological model and use a Peierls argument to provide 
heuristic basis in favour of its thermal stability. In Q 
this was proved rigorously using the formalism of quan- 
tum dynamical semigroups. In this paper we generalize 
the latter proof. We give stability conditions that apply 
to general stabilizer codes, and then apply them to toric 
and color codes. It should be noted that in the mean- 
time Bacon Q constructed a 3D model which may be 
thermally stable, too. However, our tools cannot be ap- 
plied, because the Hamiltonian is not diagonalizable in 
that case. In an independent development, the problem 
of stability of Bacon model was recently considered . 

Most likely there is no stable quantum memory in 2D. 
A full proof has been given in Q in the case of 2d Kitaev 
model (see also [11] in this context) . The no-go results for 
stabilizer codes have been provided in [1, 0] . Even though 
these results do not lead to a full proof of instability in 
the weak coupling limit, they point out quite clearly that 
there is no hope for stability in such models. Also, the 
2D Bacon model was shown instable in Q. 

Finally, let us note that many works have studied the 
existence of topological order in finite temperature (see 
e.g. [13, 0, Il9| ) and the main conclusions were compat- 
ible with the results on stability of quantum memories. 



III. SUBSYSTEM APPROACH TO 
FAULT-TOLERANT QUANTUM COMPUTING. 

The most basic approach to quantum error correction 
is the following. One considers a subspace, the code, 
chosen in such a way that, if the initial state belongs to 
the subspace, and if not too many errors hit the system, 
there exists a quantum operation, the recovery map, that 
restores the initial state. 

Such a scheme assumes ideal encoding and decoding, 
but in reality all operations are exposed to decoherence. 
The prepared state is already erroneous to a certain ex- 
tent. Then, in fault-tolerant computing schemes, it is 
subjected to repetitive error correction that keeps the 
amount of errors stable. At the end of computation, if 
the ideal recovery existed, one could map the state back 



to the code subspace, but this is not really available nor 
needed. Rather, one wants to measure some observable 
on the state. The measurement is nonideal, but if the 
state did not suffer too many errors and the measure- 
ment does not produce too many either, then the statis- 
tics will be arbitrarily close to those expected in the ideal 
scenario, where encoding and recovery are perfect oper- 
ations. 

To describe the above situation it is more convenient 
to consider a virtual subsystem (see also [2l[ in this 
context). Then the fact that we cannot apply perfect op- 
erations does not exclude that we can prepare the sub- 
system in a given state with arbitrarily high fidelity. For 
example, in standard hard drives, the sign of any single 
spin varies in time. However, the majority of signs of all 
spins (i.e. magnetization) is preserved with tremendous 
perfection. 

The virtual subsystem is actually strictly related to the 
recovery procedure (see |22| in this context). Namely, the 
observables acting on the code are lifted to the full space 
by subjecting them to an operation dual to recovery map 
[31 (we shall present this construction in section HV Cp . 

In the subsystem picture, if we were able to start with 
the state within the code the initial state would be ipv ® 
02i„c! where V denotes the virtual subsystem and the 
second subsystem is Hauc, where ^p is the wanted state 
of the subsystem. Now, the errors (if not too many of 
them occur) affect only the other subsystem, which turns 
into a mixed state. 

In reality the encoding is not error-free, so that the 
initial state is of the form \tp){4'\v^P°Anc- Then, in fault- 
tolerant schemes, there are two contradictory forces: the 
environment adds errors (which propagate), while the 
user tries to correct errors. When these forces are bal- 
anced we get a dynamical equilibrium: the system is al- 
ways in a state of the form !■(/;) (-01^ ig) PAnc- 

The above formulation is natural not only for tradi- 
tional fault-tolerant error correction schemes, but also 
for self-correcting models. In these models there is a self- 
Hamiltonian that causes the environment to do both jobs: 
adding and removing errors (the latter due to Boltzmann 
factor). In such situation, the system is (after initial 
equilibration) for a long time in the state ® Pi3i 

where pp is thermal equilibrium state of the the second 
subsystem. The Gibbs state of the whole system would 
be proportional to ly ® Pfi, and the whole system will 
eventually relax to this state. However, it may exist a 
critical temperature below which the state of V will be 
unaffected for a time that scales exponentially with the 
size of the total system [2| • 

In standard fault-tolerant schemes (where there is ac- 
tive error correction), entropy is produced by the envi- 
ronment and reduced by error correction. In the self- 
correcting case both are done by the environment. The 
fight (as usual in phase transition phenomena) is between 
entropy and energy - entropy causes errors and energy 
forces correction via Boltzmann factors. 

Let us emphasize, that we do not deal here with the so 
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called 'noiseless subsystem" [2^. The noiseless subsys- 
tem is protected against some class of errors for arbitrary 
state of the rest of the system, while here we have protec- 
tion only when the rest is in some particular state, e.g. 
Gibbs state, in the case of self-correcting models. Also 
an error correcting code usually does not form a noiseless 
subsystem: assuming that the the code can correct up to 
k errors, only if wc start with initial pure state iliv®4>Anc 
the subsystem can indeed tolerate k errors. Indeed to tol- 
erate these errors we need to start from the code, which 
in subsystem picture is the span of '0y ® ipAnc for all ip 
from subsystem, but for a fixed ipAnc- 

The sketched subsystem picture is especially useful if 
we want to analyse the interaction of the quantum mem- 
ory with the environment. Then there are known tools to 
work with autocorrelation functions of observables. To 
estimate stability of the memory it sufhccs to show that 
these autocorrelations do not decay in time. The suit- 
able observables are those from the algebra of the virtual 
subsystem. 

In section IIVI we shall describe, first in the standard 
picture, the error correction based on stabilizer formal- 
ism. Then we shall construct suitable 'dressed' observ- 
ables which define the virtual subsystem. 

IV. STABILIZER CODES AND TRANSVERSAL 
GATES 

A. Stabilizer codes 

An error correcting code [H, [l^ is a subspacc C of the 
Hilbcrt space Ti. that represents a quantum system. The 
code subspace C is such that any quantum information 
stored on it can be recovered after it has suffered errors 
from a set £ of correctable errors. Usually error correct- 
ing codes are subpaces of systems of n qubits, that is, 
Ti. = 7if " with H2 a two-dimensional Hilbert space with 
orthonormal basis {|0), |1)}. If C encodes k qubits, that 
is, has dimension 2*^, and corrects any error with support 
on less than d/2 qubits, wc say that the code is a [[n, k, d]] 
code. The support of an operator is composed of those 
qubits in which it acts nontrivially. In this context the 
size of the support of an operator O is usually called its 
weight and denoted \0\. 

An important class of error correcting codes is that of 
stabilizer codes [25l. [26|. which are defined in terms of an 
abelian subgroup S d V of the Pauli group P of n qubits. 
This is the group 

P= (i/,Xi,Zi,...,X„,Z„), (1) 

which has as generators the usual Pauli operators Xi, Zi 
on the i-th qubit, 

^ = l+>(+|-|->(-|, Z=|0)(0|-|1)(1|, (2) 

where |±) -.^ (|0) ± \ l))/V2. The subgroups S is called 
the stabilizer of the code C because its elements £ C 



arc defined by 

Vs€5 s|V') = |V'>- (3) 

Note that S cannot contain the clement — /. The number 
of encoded qubits in a stabilizer code is 

k^n-g, (4) 

with g the number of independent generators of S. To 
give an expression for the distance of the code we have 
to introduce the normalizcr group Af. This is the nor- 
malizcr of S in V, which coincides with the ccntralizer 
of 5* in V: its elements n G J\f C V commute with all 
the stabilizers s G S. Then the distance of the code 
C is given by the minimum weight of the elements of 
Af — S, which indeed implement nontrivial unitary op- 
erations on encoded states. In fact, it is possible to 
choose among them the X and Z Pauli operators for 
the k encoded qubits: Xi, Zi, . . . , X^, Z^. € A/" — 5. 
These operators must satisfy the usual Pauli algebra, 
so that we_have Xf = Z} = 1, {Xi,Zi\ = and 
[X„ Zj] = [X„ Xj] = Zj] = for ^ ^ j. 

B. Error correction 

If we put some quantum information into a given code 
subspace we know that as long as the system does not 
suffer too many errors it is in principle possible to recover 
the information without any losses. But of course we 
are interested in recovering it in practice, which is not 
necessarily just as easy. 

Stabilizer codes have a particularly simple structure 
that makes recovery from errors relatively easy. The first 
step in error correction is the measurement of a set of 
generators Si of the stabilizer 

5 =< Sl, . . . ,Sg > . (5) 

Let us denote the eigenvalue of the i-th generator after 
the measurement as (—I)''', 6^ = 0, 1 G Z2. The list of 
eigenvalues b = {hi) E is called the error syndrome. 
The second step is then to choose, using the syndrome s, 
a Pauli operator K G V such that 

Ks, = (6) 

We set synd(iir) := b when ([6]) holds, so that synd defines 
a group homomorphism 

synd : V — > Zf (7) 
K^h, (8) 

with kernel Af. 

If the correction operator K is applied to the system 
it will return to the code subspace C. Whether it returns 
to the original state that we had encoded will depend 
on the errors that occurred and the choice of K. To 
clarify this point, suppose for simplicity that the action of 
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errors amounts just to the application of a Pauli operator 
E G V to the system. Then we have by assumption 
KE e M, and the error correction procedure will succeed 
if in addition KE G S. 

In summary, the whole error correction procedure is 
encoded in a function 

corr : Zf — yV (9) 
b I — y K, (10) 

that takes syndromes b to correction operators K . This 
corr is such that 

synd(corr(b)) = b. 

One should choose corr in such a way that the probabil- 
ity of success for a given source of errors is maximized. 
However, in practice it is also important that corr can 
be computed efficiently (say, in polynomial time in the 
number n of physical qubits), a requirement that may 
reduce the success probability with respect to the ideal 
case. 



C. Dressed observables and critical syndromes 

Take any nontrivial self-adjoint encoded operator N G 
Af — S. Suppose that we initialize the system in an en- 
coded state and we want to keep track of the value of 
these encoded operator as errors occur, from a theoreti- 
cal perspective. In particular, we would like to keep track 
of the value of N were we to correct the system using the 
above procedure. We may for this purpose introduce the 
operator 

TVd, ^ corr(b)^7Vcorr(b)Pb, (H) 
bezf 

where Pb is the projector onto the subspacc with syn- 
drome b, that is. 

9 
i=l 

Note that for encoded states iVdr = N. The operator 
A^dr is a more stable version of N, because as long as 
the system only suffers correctable errors its expectation 
value does not change. This is not the case for A'', which 
may change even with an error on a single qubit. 

Given an encoded operator N it is useful to introduce 
the notion of critical syndromes as those which are just a 
single-qubit error apart from changing the value of A^dr- 
Let \ipb) denote a normalized state with PblV^b) = |^b) 
(there are a few of them, but what follows does not de- 
pend on which one we will choose). 

For any Pauli operator E G V we have 

iVdrPl^Z-b) ^iSiVdrlV-b). (13) 



In particular, the sign is negative if and only if the en- 
coded operator 

corr(V'b) E corr(V'b + synd(i?)) (14) 

anticommutcs with A^. For later convenience we set 

SiN, E,h):^ {iPk I nI E^ A^dr ^ | V-b) ■ (15) 

Notice that for E = E1E2 we have 

S{N,E,h) = SiN,Ei,h + synd{E2))S{N,E2,h). (16) 

When there exist a singlc-qubit Pauli operator a such 
that S{N,a,h) ~ —1, we say that b is A^-critical. For 
each A^, we denote by cntN C the set of A^-critical 
syndromes b. Similarly, we denote by crit'^ C P x Zj 
the set of pairs {E,h) with ^(A^,^;, b) = -1. Due to 
P^ . if {E,h) g crit^ there exists a decomposition E = 
E1E2 such that the Ei-s have disjoint supports and b + 
synd(i?2) G critAr. 

D. Transversal gates 

An important feature of several classes of stabilizers 
codes is that they allow the implementation of logical 
gates transversally, that is, performing unitary gates on 
individual qubits or in small disjoint sets of them. Usu- 
ally we are interested in stabilizer codes C that encode a 
single qubit. Then a one-qubit logical gate G that can 
be implemented in C through a unitary operation of the 
form 

U = ^U, (17) 

i 

is transversal. Here the index i runs over physical qubits. 
Frequently one has [/; = Uj , which means that no specific 
addressing of the physical qubits is necessary, they are all 
treated on equal footing. 

A particularly simple case is that of U belonging to the 
normalizer of V, usually called the Clifford group. Then 
U is a. tranversal gate if for any s e 5 

UlsU e S, (18) 

because for any G C we must have sU\^p) = U\ip). An 
advantage of such tranversal gates is that, when applied 
with perfect accuracy, they affect the error syndrome of 
the code in a definite way. More generally it is enough 
for U to satisfy 

U^PoU = Pq. (19) 

In section IVirF3l we will consider such gates, which are 
essential for universal quantum computation. It is clear 
that from condition (|19[) we get no information about 
how U affects the syndrome. However, (|17p is a big con- 
straint in this regard. In particular, errors cannot spread. 
That is, if is an encoded state and E an error, then 
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UE\ip) = E'U\ip) for some E' with the same support as 
E. A less trivial observation is that, in a certain sense, 
negative error syndromes also do not spread. Suppose 
that we choose a possibly overcomplete set of stabilizer 
generators which are local in a given lattice of a certain 
dimension where the qubits live. Take a region R of this 
lattice such that no encoded operator has a support com- 
pletely contained on it. Then if O is a Pauli operator that 
commutes with all the local generators that have part of 
their support in R, it follows as a slight generalization of 
the cleaning lemma of that 0\i;) = O'lV') for some O' 
with support only outside R. This means that, if E\ip) 
is free of negative error syndromes in a region, then the 
same will be true for UE\tp) because i?|V') = E"]!/;) for 
some E" with support outside R. 

Similarly we can consider transversal two-qubit logical 
gates, which are typically implemented on a pair of equal 
stabilizer codes through unitaries of the form p?)) . where 
now i runs over pairs of equivalent physical qubits on the 
two codes, so that the Ui are two-qubit gates. Finally, 
it is also possible to consider transversal measurements, 
where each physical qubit is measured on a given basis 
and from the net result a measurement for the encoded 
qubit is obtained. 

E. CSS-like codes 

CSS-like codes are those stabilizer codes for which the 
generators of the stabilizer can be chosen so that S = 
(sf ,...,s^,sf,...,4), 31 +.92 =g, 

sf e Vx, sf e Vz, 

with Vx (^1, ...,X„),Vz ■■= {Zu ■■■,Zn). In such 
codes, X-typc (Z-type) logical operators may be chosen 
from Vx (J-'z), that is, so that 

X, e Vx. Z, e Vz- 

Then X-typc (Z-type) errors can only change the value of 
Z-tyj>e (X-type) encoded operators, and their correction 
only involves Z-type (X-type) stabilizers. That is, the 
error correction procedure can be divided into two, and 
the relevant functions are 

synd^ : Vz Zf \ corr^ : Zf^ Vz, 
synd^ : Vx Zf , corrx : Zf Vx- (20) 

Then for any Z-type encoded operator N C JV O Vz the 
corresponding dressed operator takes the form 

A^dr:= J2 coTVxihy NcoTTx{h)P^, (21) 
hezp 

where is the projector onto the subspace with Z- 
syndrome b, that is, 

92 

■■=2-^^-11 {I + {-If ^sf). (22) 

i=l 



The dressed X-type Pauli observables are analogous. 

An important consequence of the form of the dressed 
observables pT|) is that transversal measurements are 
possible. Namely, if we measure all physical qubits in 
the X basis {Z basis) we can recover the value of all 
encoded X-type (Z-type) operators because error cor- 
rection commutes with measurements. Indeed, after the 
measurement correction is entirely classical. 

Finally, an analogous transversal initialization is pos- 
sible. If we initialize a product state of the form |0)®" 
(|-l-)®") and measure all X-type (Z-type) stabiHzer op- 
erators, the resulting state has well defined Z-type (X- 
type) encoded operators. Namely, Z^ = 1 (respectively, 
Xi = 1). The state may have many X-type (Z-type) sta- 
bilizers with the wrong sign, but indeed the correction of 
Z-type (X-type) errors is immaterial because such er- 
rors cannot change the encoded state. Of course this is 
the case when environment is absent. We shall discuss 
initialization in the presence of noise in section IVIII Al 



V. MODELING INTERACTION WITH 
ENVIRONMENT 

In this section we shall describe the evolution of an 
open system in the weak coupling approximation. Next 
we shall recall the relation between the fidelity criterion 
and the decay rate of observables, and a useful upper 
bound for the latter Q (see also [23] in this context). 



A. Evolution in weak coupling limit 

We consider a quantum system with discrete energy 
spectrum is coupled to a collection of heat baths leading 
to the global Hamiltonian 

i/ = iT^y ^ + '"^^'H 7? with if = ^ S'a ® /a , 

a 

(23) 

where the Sa are system operators and the fa bath op- 
erators. 

We assume that the baths are independent, which ex- 
cludes such phenomena as decoherence free subspaces 
(see e.g. Hfl [M Hi]). The evolution of the system in 
Heisenberg picture in weak coupling limit is the follow- 
ing 

dX 

— =i[H'y',X]+CdUX)=:C{X) (24) 
at 

>Cdis(X) - ^ E E ha{^){siH [X, Sa{Lo)]+ (25) 
+ [SiicolX]Saico)) (26) 

Here ha are Fourier transforms of the autocorrelation 
functions of the fa- They describe the rate at which the 
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coupling is able to transfer an energy lo from the bath to 
the system. In the case of thermal baths, they satisfy 

K{-uj)^(r^^K{uj) {27) 

which is a consequence of the KMS condition [s^l- The 
dissipative part commutes with Hamiltonian part, i.e. 
[H,C] = (where H{X) = Therefore in 

analysis of thermal stability, it is enough to consider 
solely the dissipative part. An important property of 
the dissipative generator is that — £ is positive definite 
in scalar product 

{X,Y)0:=TyppX^Y. (28) 

Moreover we have £(1) = 0, and if the commutant 
of the system operators Sa and H'''^^ is trivial, then the 
eigenvalue is nondegenerate. This means, that there are 
no constants of motion - all other observables eventually 
decay to identity. 

B. Fidelity and autocorrelation functions 

Since the generator is hermitian in scalar product ()28p . 

a strong tool is spectral analysis of C. Namely, for any 

observable satisfying {A^A)p = 1 and {A,I)p — one 
easily finds that 

{A,A{t))p>e-^' (29) 

provided that 

-{A,C{A))p<e. (30) 

Thus, to show that an observable A is stable, it is enough 
to estimate —{A,C{A)), which can be therefore called 
decay rate for the observable A. If this quantity decreases 
exponentially with size of the system, we obtain stability. 

However, the criterion of stability of quantum infor- 
mation is measured in a most direct way by fidelity. Yet, 
in some natural situations, the fidelity criterion is equiv- 
alent to stability of two complementary observables in 
terms of autocorrelation function. Namely suppose that 
we divide the system into two: Tiy ® Ti-Anc, where the 
first system is the qubit to be protected. Let X and Z 
be two Pauli observables, which act on the memory, i.e. 
they are of the form Xy (8) Iauc and Zy CS" Iauc- Suppose 
further that the Gibbs state is proportional to identity on 
the system V we have [H the following bound for average 
fidelity 

F > i((X, e^'X)p + {Z, e'''Z)p) > e"^* (31) 

where e is upper bound for the rates —{X,C{X))p and 
~'{Z,C{Z))p. Therefore, to prove stability of quantum 
memory it is enough to bound the above decay rates. 
One can actually show, that the first inequality holds for 
whatever physical operation (trace preserving completely 
positive map A, i.e. we have 

F>^{{X,AiX)) + {Z,AiZ))). (32) 



C. Upper bound for decay rate 

In Q the following useful bound for decay rate was 
obtained: 

^{A,C{A))p < 2/l„ax^([^a,A],[^a,A])^ (33) 
a 

where 

ftmax = sup ha{uj). (34) 

a,Lj>0 

The advantage of this estimate is that it refers directly 
to the system operators Sa rather to their Fourier trans- 
forms. 



VI. CONDITIONS FOR SELF-CORRECTING 
QUANTUM MEMORIES 

A. Stabilizer code models 

Given a stabilizer code, we may consider a Hamiltonian 
model of the form 

r 

H = -Y^t,s[, (35) 

where ti > are coupling constants and {s^} is a possibly 
overcomplete set of r generators of the stabilizer S. The 
elements of the algebra of encoded Pauli operators gen- 
erated by Xi, Zi, . . . , Xk, Zk commute with H, showing 
that each energy level is at least 2'^-fold degenerate. For 
a particular choice of the encoded operators, one may 
split the Hilbert space 

n = Henc ® T^oxc = C^'^ ® C^f , (36) 

in such a way that Ticnc represents the encoded degrees 
of freedom and Tioxc characterizes excitations. Let us 
emphasize that the above division is not yet the one we 
considered in section IIIII Indeed, the subsystem Ticnc is 
defined by the bare observables, and it is not protected. 
In particular, we may choose a basis |b) in Tioxc with 
the elements labeled by b e so that setting \ip,h) := 
I?/;) (X) |b) for any £ TYcnc we have 

5,|V;,b) = (-l)''-|V^,b). (37) 

Recall that Si are generators from the complete set ([5|). 
Then the states of the form \tp,h) are eigenvectors of the 
Hamiltonian (j35[) . namely 

r 

if IV, b) = - ^ tf b) =: Ey,\^, b), (38) 

where Ci(b) = 0, 1 are obtained from 
g 

^^^n^f' c,,=0,l, (39) 
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as Ci(b) = ^ACijhj, with addition modulo 2 



Gibbs state for the Hamiltonian (133 



Thus the 
can be written as 

(40) 



we will derive for a wide class of systems a set of suffi- 
cient conditions for this stability to happen at some finite 
temperature. 



B. A bound in terms of critical syndromes 

In this section we shall show that the decay rate for a 
dressed observable iVdr is bounded by the probability of 
a critical syndrome multiplied by some polynomial in the 
size of the system. Thus, if the probability of the critical 
syndrome is exponentially small in the size of the system, 
the observable is stable. 

Suppose that we have a Hamiltonian model psp that 
interacts with a bath through a coupling of the form 

H^nt^Y^X,® f, + ^Y,®g, + ^Z,®h^. (41) 

i i i 

For simplicity, we will write 



int 



(42) 



with (T running over the three types of Pauli operators 
and over all sites. We will suppose that our model enjoys 
the following two properties. First, the ti couplings take 
only a finite number of values. Let tmax be the maxi- 
mum of these values and imin the minimum. Second, the 
number of stabilizers Si (from psp ) that contain in their 
support a given qubit is bounded. When these conditions 
hold, it follows that the number of different lj-s is finite 
and independent of the size of the system. Therefore 
/imax is finite and (|33p applies. 

Let us go back to the dressed operators introduced in 
section HV CI It is our aim to check how stable they are 
under thermal fluctuations. From ([55)1 and we have 
that for any encoded operator N the decay rate of iVdr 
is bounded as follows 



{N,\,,C{N^,))l3 < 



< 8K 



El 



E,(l- W^,b)) 



< 



(43) 



Here n is the total number of sites, which we will suppose 
that grows polynomially with the size of the system L, 
and Pcritjv is the total probability at a given temperature 
of finding a syndrome that is critical for N, that is 



E 



critjv 



bGcriti\ 



El 



-/3Sb 



(44) 



When for a given temperature Pcritjv is exponentially 
small on the size of the system L we say that TVdr is 
stable: it is self-protected by a free energy barrier pro- 
duced by the Hamiltonian ([55]) . In the following section 



C. A criterion based in connectedness 

Here we present a set of conditions which assures that 
the probability of critical syndrome is exponentially small 
in the size of the system. These conditions require the ex- 
citations to be connected in a suitable way forming clus- 
ters, and also that critical configurations always contain 
clusters of size comparable to the system size. Thus it is 
essentially an application of Peierls argument [lU, [H, HI] 
as first done in j2| in the context of quantum computa- 
tion. 

Suppose that the set {s^} of generators of the stabi- 
lizer iS that appear in Hamiltonian (j35[) are the nodes 
of a graph F. Alternatively, the nodes represent exci- 
tations. Then any syndrome b g corresponds to a 
subset nds(b). Namely, the i-th node belongs to the set 
nds(b) if Cj(b) = 1. We may introduce paths in the graph 
in the usual way. This gives rise to a notion of connect- 
edness: a set of nodes is connected if between any two of 
its nodes there exists a path such that all its nodes are 
contained in the set. We say that b and b' are connected 
if nds(b) and nds(b') are not disjoint or if any of the ele- 
ments of nds(b) is linked to an element of nds(b'). Also, 
we say that a site j (or physical qubit) touches a node s[ 
in F when either {X^, s-} = or {Zj, s-} = 0. 



Rather than a single Hamiltonian model psp . we con- 
sider a family of models with increasing system size 
L, each with a corresponding graph F^. We assume 
that |Fl|, the number of nodes, is polynomial in L. The 
graphs Tl must satisfy the following conditions for cer- 
tain integers A, fi, v and a positive number ^ (all indepen- 
dent of L): 

(i) The number of nodes touched by a site is at most A. 

(ii) The number of links meeting at a node is at most \x. 

(iii) For every \i ^7j\ there exists a decomposition 



(45) 



such that the b^-s are not mutually connected and have 
at most V connected components each, 
(iv) If |nds(b)| < then b ^ critw for any encoded 
operator TV and given any b' not connected to b 



corr(b + b') = corr(b) corr(b'). 



(46) 



The goal is to give a bound for the probability Pcritw ■ 
We first want to bound the probability Pbo to find an 
excitation configuration nds(b) such that bp is one of 
the components in the decomposition (j45|) . Denoting by 
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{b'j,} the set of all such configurations we have 



b„ 



GZ; 



3 e 



s e 



< 



We now show that, for a given b and a encoded op- 
erator N, h £ critM only if it has a component bg 
with |nds(bo)| >= ^L/X^. Suppose to the contrary that 
b = El=i t»i with |nds(bi)| < ^L/Xfi. If tr is a single- 
qubit Pauli operator, then (i,ii) imply that at most Xfi of 
the bi-s are connected to nds(synd(CT)), say bi, . . . , b^. 
Using (iv) we get 

r 

corr(b) = corr(b') corr(b,;), 

r 

corr(b + synd(cr)) = corr(b' + synd((T)) corr(bi) 

i=fc+l 

(48) 

with b' = bj- It follows from (HH) that S{N, cr, b) = 

S{N,a,h') and thus (iv) that b ^ critjv. 

Finally, it is easy to show that the number N{1) of 
connected collections of nodes of size I that contain a 
given node is bounded as follows (see appendix [X)) : 



Nil) < e 



tI 



(49) 



where t = yLtln2. Thus, the total number of collections of 
I nodes with at most v connected components is bounded 
by iFLl'^r-ie^' < jFLp'^e^'. 

Putting together these results we get 



1 



1 - e- 



Pcrit„ < iPip'' e-*' = poly(L)e 

(50) 

with 6 = /3tmin — T positive below a critical temperature. 
This shows the stability of N^t under the given assump- 
tions. 



D. Resilience to noise 

It is natural to expect that a self-protecting quantum 
memory will also be resilient, to some extent, to sud- 
denly applied noise. This is interesting since it allows us 
to switch off the protecting Hamiltonian for short times, 
since the unprotected period is equivalent to a noisy pro- 
cess. As we shall show, as long as the amount of noise 
is below a threshold, we are safe. This threshold will 
reduce as we increase the temperature till it vanishes at 
the critical temperature. 

To fix ideas, let us consider a process Tp in which each 
qubit undergoes a depolarizing channel characterized by 



an 'error' probability p. That is, any given operator TV 
is mapped in the Hcisemberg picture to 



Tp(iV) = ^(1 - p)"-l^lpl^l E'^NE, 



(51) 



where the sum runs over a set of representants of V jil. 
We arc under the error threshold when in the L oo 
limit for any given encoded operator N we have 



But, setting rj 

(iVdr,Tp(7Vd..))/3 



log((l-p)/p), 

Ebije-'3''''""""^(^dni?,b) 



That is, {Ndr,TpiNd,))p 



Eb.£;e-'^^-"l^ 
= 1 



(52) 



(53) 



2-Pcrit' with 



crit' 



E(b,B)ecrit^ 



-/3Et,-ri\E\ 



Eb.ije-^^"-"!^! 



(54) 



This expression is similar to p4|) , and we can use a sim- 
ilar reasoning to show that Pciit^ ^ for low enough 
temperature and noise. 

First, we attach to each Pauli operator E the set 
nds(£^) C Fi consisting of all nodes touched by any site 
in the support of E. Notice that nds(synd(ii^)) C nds(£') 
and |nds(i;)| < A|£;|. Thus, when A|£;| + |nds(b)| < S,L 
we have (£',b) ^ crit'^. 

Secondly, we need a substitute of the decomposition 
PS)) to deal with pairs {E,h). To this end, we set 
E = Y[j Ej in such a way that the sets nds{Ej) form the 
connected components of nds(£'). Then we arrange the 
bi-s and Ej-s in maximal connected collections. That is, 
if one of these collections is {b^j^ , . . . ba^ , Ea^ , . . ■ , Ea^ }, 
then the set nds(bai) U • • • U nds(baj U nds(£:ai) U • • • U 
nds{Ea^) is connected and disconnected from any set 
nds(bi) or nds(Ej) that is not part of the collection. We 
denote each collection as {Ea,ha), with ha = Ei'^iiJ 
Ea = E^ Ea, , SO that b = and E = Y{^Ea. This 

gives the desired decomposition of {E, b) . 

Finally, we observe that when X\Ea\ + |nds(ba)| < 
for all pairs (i?a, ba) we have 

corr(b) = JJcorr(ba), 

a 

corr(b + synd(^;)) = ]J corr(ba + synd(^;a)), (55) 



so that S{N, E, b) = Ha ^i^, Ea, ha). We omit the rest 
of details. 



VII. TOPOLOGICAL STABILIZER CODES 

In this section we shall construct codes that allow a 
universal set of transversal gates and, at the same time, 
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give rise to self-correcting local Hamiltonian models. The 
codes will belong to the important class of local stabilizer 
codes, in particular topological codes. These are con- 
structed from manifolds where the physical qubits are 
placed, and have two main characteristics: 

• There exists a set of generators of the stabilizer S 
with a spatially local support. 

• The encoded operators in Af — S are global, they 
have a support which is topologically nontrivial. 

We will give below a more precise characterization. The 
quantum Hamiltonian models (j35|) obtained from topo- 
logical codes are topologically ordered: their ground 
spaces have a degeneracy of topological origin and their 
excitations are gapped and subject only to topological 
interactions. As we will see, excitations can be either 
particles or higher dimensional objects. Self-correction 
will appear in codes with excitations of the latter kind. 

The original example of topological stabilizer codes 
are toric codes, as introduced by Kitaev[l|, These 
codes and their generalizations arc based on homology 
theory [l^l- Indeed, the number of encoded qubits in a 
given manifold is dictated by its Betti numbers, i.e., the 
dimensions of homology groups, which are topological in- 
variants. 

An important drawback of toric codes is that they do 
not allow the transversal implementation of a sufhciently 
rich set of gates. There exists however a class of topolog- 
ical stabilizer codes that surmount this difficulty: color 
codes. In 2D Q they allow the transversal implemen- 
tation of the so-called Clifford group, the normalizer of 
the Pauli grou p in the group of unitary operations on n 
qubits. In 3D [lO| they allow the transversal implemen- 
tation of a set of operations that, when measurements 
are included, gives rise to transversal universal quan- 
tum computation. Below we will analyze in detail the 
transversality properties of general £)-dimensional color 
codes. 



A. Local stabilizer codes 

A family of codes {Cl} is local if there exists integers 
m and n such that 

(i) each stabilizer Sl has a set of generators {si, . . . , Sg^^} 
such that \sj\ < a, j = 1, . . . , (7l, 

(ii) each physical qubit is in the support of at most /3 
such generators Sj and 

(iii) for any integer d, there exist a code Cl with distance 
bigger than d. 

Given a code Cl in the family, it is always possible to 
construct a local quantum Hamiltonian model such that 
its ground state is Cl- Namely, we may set 

Hl = -J2^^- (56) 

4=1 



Topological stabilizer codes are local codes in which 
locality has a geometrical meaning. Physical qubits are 
placed in a lattice in a given manifold, in such a way that 
the generators of the stabilizer have a support contained 
in a small neighborhood of the lattice. The support of 
nontrivial encoded operators, in turn, must be global, in 
such a way that the distance of the codes grows with the 
size of the lattice [s^ . 

Below we will study in full generality the two main 
examples of topological stabilizer codes, toric codes and 
color codes. We will unify the treatment of both types of 
codes through homology theory. 

B. Homology groups 

Recall that a homology group is defined from an exact 
sequence such as 

A,^A^^ Ao (57) 

where Ai are Abelian groups and 9i, the boundary oper- 
ators, are homomorphisms with 

910^2 = 0. (58) 

The elements of the subgroup Zi := kerc^i C Ai are 
called cycles. The elements of the subgroup Bi 
im92 C Ai are called boundaries. Then the condition 
([55)1 states that all boundaries are cycles. Bi C Zi. Be- 
cause of this, we can divide cycles in homology classes. 
Two cycles z, z' are homologous, z z' , if z — z' is a 
boundary. The homology group H has as its elements 
the corresponding homology classes. It is the quotient 

H:=Zi/Bi. (59) 

In our case, the elements of the groups Ai will be chains 
of suitable geometrical objects from a lattice in a given 
manifold, and the homology group will be a topological 
invariant of the manifold, independent of the particular 
lattice considered. We are only interested in Z2 homol- 
ogy, in which the chains in Ai arc binary formal sums of 
geometrical objects in a set Si. That is, if a e Ai then 

a=^6,s, 6, = 0,1. (60) 

When needed, we will regard Z2 chains as sets in the 
usual way. That is, we identify the chain a above with 
the set {s e |6s = 1}. 

When working with such Z2 homology groups we may 
consider the dual exact sequence 

Ao ^ Ai ^ A2. (61) 

The dual boundary operators d* are defined so that for 
a e Si-i, b G Si we have a G dib if and only if 5 e 
i9*a. The dual exact sequence (|5T|) gives rise to a dual 
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homology group H* := Z^/Bl which is isomorphic to H. 
Notice that for a £ ^i, S2 € 5*2 and Sq G Sq we have 

So S dia <;=4> \dlsQ n a| = 1 mod 2, 

S2 e d^a 19252 n a| = 1 mod 2. (62) 

C. Codes from homology groups 

Given a Z2 homology group H obtained from a lattice 
in a manifold, as those discussed in the previous section, 
we can introduce in a natural way a CSS-like stabilizer 
code C with a basis labeled by the elements of H. The 
procedure can be summarized as follows. Recall that the 
abelian groups At in the exact sequence (f57|) are con- 
structed from sets of geometrical objects Si. To define 
the code, we attach physical qubits to the elements of 6*1 , 
X-type stabilizers to the elements of S2 and Z-type sta- 
bilizers to the elements of 5*0. In what follows, we detail 
the procedure. 

First, we attach a physical qubit to each of the elements 
of 5*1 = {rj}, so that n = \Si\. For each a G Ai with 

a = Uj Vj, aj = 0, 1 (63) 
j 

we define a X-type and a Z-type Pauli operator setting 

Xa:^(^X;\ Z,:==(g)Z;^ a = ^c,r^. (64) 
j 3 

Note that the definitions (|64p describe two group isomor- 
phisms. Moreover, 

XaZa' = {-lf'''\Za'Xa. (65) 

Then using ^ 

a e Zi ^ V6 e B*, [Xa, Zb] = (66) 
a e Z* ^ V5 e Bi, [Za, Xb] = 0. (67) 

As a consequence of ((55)) . [Xf,, Z^'] = and wc can define 
the stabilizer group S as the subgroup of V with elements 
of the form 

XbZb', beBi,b'eBl (68) 

It follows then that the operators of the form 

X,Z,,, z£Zi,z'eZ*, (69) 

form, up to phases ±l,±i, the normalizer group Af. If 
Xi, . . . ,Xk G Vx and Zi, . . . ,Zk G Vz form a Pauh basis 
for encoded operators, then Xi = X^^ and Zi ~ Z^' for 
suitable cycles g Zi, z- g Z*. The cycles Zi (z-) are 
representants of classes Zi G H (z^' S H*) that form an 
independent set of generators of H (H*), and thus wc 
have the desired result. 



It is instructive to work out explicitly a basis for the 
resulting code C. To this end, for any a € Ai we set 

|a)o XJO)®", |a)+ :==ZJ+)®". (70) 

Then {|a)o} and {|a)-)-} are bases of the Hilbert space, 
with elements labeled by the chains a G Ai. For any 
cycles z G Zi, z' G ZJ' we define the states 

|z)o:=^|c)o, |z')+ := ^ |c) + . (71) 

Then, as can be easily checked, {|z)o} arid {|z') + } are 
bases for the topological stabilizer code C with elements 
labeled by homology classes of H and H* , respectively. 

1. Error correction 

We analyze now error correction in the code C just de- 
scribed. First, we need a generating set for the stabilizer 
group, which can be obtained by attaching an operator 
to each of the elements of S'2 and Sq. Set for any 02 G ^2, 

oo G Ao, 

'■= Xd2a2: ^ao '■=Zo*ao- (72) 

Then if S'2 = {pi, . . . Sq = {gi, . . . , g,„}, we may 
write 

S :~ {Xp-^ , ■ ■ ■ , Xp^ , Zgj , . . . , Zq^). (73) 

Indeed, the above set may be overcomplete. We can 
choose in that case suitable subsets S'2 = {pi , ■ . . , Pg^ } C 
S2, S[ — {qi, ...,(733} C Si. We may naturally regard 
any S2 C S2 as some b G Zf^ , setting bi ~ I when pi G S2, 
and similarly for S[ and Z2^. Then, for a G ^1, 

synd^(Za) = S[ D d^a, 

syxidz{Xa) = S'^r^dia. (74) 

That is, for a given Pauli error XaZa' the error syndrome 
gives us information about the boundaries of a and a'. 
Error correction amounts to choose an operator XcZ^' = 
corr(syiid(XaZa')) with dic = dia and 9|c' = 5|a'. It 
will succeed if and only if a+c G Bi and a' +c' G B*, that 
is, if and only if the correction operator and the error are 
equivalent up to homology. 

2. Hamiltonian and excitations 

Given the generators of the stabilizer in (|73p . we can 
consider the quantum Hamiltonian model 

/ m 

H = -Y,Xp^-Y,Zq^. (75) 

1=1 1=1 

The excitations of this model are gapped, with two units 
of energy per stabilizer violation. We say that the excita- 
tions related to Xp^ {Zq.) generators are X-type (Z-type) 
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excitations. Then, as follows from the previous analysis 
of error correction, X-type excitations are arranged in 
the form of _B* boundaries and Z-type excitations are 
arranged in the form of Bi boundaries. 

D. Generalized toric codes 

Consider a lattice in a D-manifold. We call the n- 
dimensional elements of the lattice n-cells: 0-cells are 
vertices, 1-cells are links, 2-cells are faces, and so on. For 
n = 0, D, let Sn be the set of all n-cells, and Cn the 
corresponding abelian group of chains. We in addition 
set Cd+1 = C_i = 0, with the trivial group. For 
n = 1, . . . Z? we introduce the boundary operators 

Cn ^ Cn-1 (76) 

in the usual way. That is, they arc homomorphisms such 
that if the {n — l)-cells s,; G Sn-i form the geometrical 
boundary of the n-ceW s then ds = {si}. We also define 
the homomorphisms do+i and 9o in the only possible 
way. 

We can obtain a generalized toric code for each d with 
d — . . . ,D. This is the code that corresponds to the 
exact sequence 

Cd+i ^Cd^ Ca-i, (77) 

as described in the previous section. The corresponding 
homology group is the standard d-i\\ homology group of 
the manifold, 

H, := (78) 

Thus the number of encoded qubits k is h^, the d-th 
Betti number of the manifold. For example, for d = 1 wc 
have the homology group of curves, so that the number 
of encoded qubits is equal to the number of independent 
non-boundary loops in the manifold. Notice that the 0- 
th and the D-ih code are just classical repetition codes, 
with one encoded qubit per connected component of the 
manifold. 

The dual sequence can be visualized more easily in the 
dual lattice, in which r-cells become (£)— r)-cells. Indeed, 
in the dual lattice the dual exact sequence corresponds 
to the exact sequence of the d-th homology group, where 

d:= D-d. (79) 

Thus, in a d-th color code the relevant homology groups 
are Hj and Hd ■ Cycles take the form thus of closed d- and 
d-branes, which are mapped to operators using (|64p . The 
commutation rules of closed brane operators are purely 
topological. Namely, given a closed d-brane z £ Zd and 
a closed d-brane z' € Z^, the operators and Z^i anti- 
commutc if z and z' intersect and odd number of times. 
Otherwise they commute. 



1. Branyons 

For d= 1, . . . , D — 1, the Hamiltonian model (jTS]) of 
a d-th toric code contains two kinds of brane-like exci- 
tations, which are created by open d and d-brane oper- 
ators. Z-type excitations are closed (d — l)-branes and 
X-type excitations are closed (d — l)-branes in the dual 
lattice. Indeed, both kind of excitations are not only 
closed but also boundaries. For example, for 13 = 2 and 
d = 1 the excitations are 0-branes, that is, particles. 
They can only be created or destroyed in pairs since they 
have to be boundaries of strings. These two kinds of par- 
ticles interact topologically and are thus called anyons. 
More generally, Z-type (d— l)-branes and X-type (d— 1)- 
branes interact topologically and thus deserve the name 
of branyons [lH . 

Since (n — l)-branyons appear at the boundaries of n- 
brane operators, n = d, D — d, it follows that such opera- 
tors represent branyon processes. In particular, processes 
that do not create or destroy excitations are related to 
closed branes. Then one can derive from the commuta- 
tion rules for closed brane operators the topological ±1 
phases appearing in branyon processes. 

2. Self- correcting models 

To fix ideas, let us consider a family of hypcrcubic 
latices in a D-dimensional torus. Each lattice has the 
same lenght L in the D directions. We will show that, for 
d > 1, the models with Hamiltonian ([75)1 for the d-th code 
protect Z-type encoded operators in the sense of section 
Then it follows by duality that for d < D - 1 X-type 
operators are protected, and thus that for 1 < d < Z3 — 1 
we have self-correcting quantum memories, which thus 
exist for D > 4. 

Recall that we have to define a connection scheme for 
excitations and a correction procedure such that the con- 
ditions (i)-(iii) of section fVI CI arc satisified. The notion 
of connectedness of excitations, that is, of (d — l)-cells, 
is naturally given by the lattice. To translate this into 
the language of section IVI CI we just construct a graph 
with (d — l)-cells as its nodes in which nodes are con- 
nected by a link when the corresponding cells meet at a 
(d — 2)-cell (notice that d > 2). Condition (i) is clearly 
satisfied. As for condition (ii), we first notice that any 
connected set of excitations is closed. Some connected 
sets are boundaries, so that they can be switched off in- 
dependently and give no further problems. But those 
which have non-trivial homology cannot be switched in- 
dependently. Instead, we have to gather several of them 
together to get a boundary, but only in a finite amount 
because Hd is finite. Indeed, ?/ = \Hd\ is a trivial bound 
and condition ii is satisfied. 

Finally, we need a correction procedure, that is, the 
function corrx in (PO)) . We define it in two steps. In the 
first step we divide the syndrome b in the components 
hi of In the second step, for each we find a suit- 
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able collection c' of d-cells such that dc' = c = nds(bi). 
The choice of the collection of d-cells c' must be such 
that it ensures stability. In fact, we only consider those 
connected sets of {d— l)-cells that can be put inside a hy- 
percube of length I with I < L, and choose any c' that is 
contained in the hypercube. This task that is always pos- 
sible because c has trivial homology, and it can be easily 
checked that it can be accomplished in polynomial time 
in I. We leave the correction of other sets of (d — l)-cells 
undefined, because it is not important. 

Now that we have a correction procedure, we can check 
condition iii. It suffices to show that no Z-syndrome 
b e Z^^ that only contains connected sets of size at most 
L/2 — 1 can be critical. To this end, take any d-cell t 
and let us check whether b is Xj-critical. Let Si be the 
connected components of the set of {d — l)-cells related 
to b with which t is in contact. Both t and these sets 
Si arc contained in a hypercube of length L — 1, because 
each Si is contained in a hypercube of length L/2 — 1. 
Therefore, using the notation of ([7^ we have that the 
operator corrx(b)corrx(b + 5*2 n dit) has support in the 
hypercube of length i — 1 and thus has trivial homology. 
This implies condition iii and thus the the stability of 
iVdr for any given N C Vz- 



E. Generalized color codes 

Color codes are local codes that are constructed in a 
special kind of lattices called D-colexes that live in D- 
manifolds. We will first describe color codes without re- 
ferring to homology, and then give a homological picture 
to fit them in the general framework given above. As in 
the case of toric codes, several color codes labeled with 
d = 0, . . . , Z? can be obtained from a given lattice, which 
again will display d—l dimensional and D — d—1 dimen- 
sional excitactions. However, unlike in toric codes these 
excitation can branch and form nets. 



1. D-colexes 

A D-colex A is a _D-dimensional lattice with certain 
colorability properties. Such lattices A can be described 
(m in terms of their 1-skeleton or graph, wich is 1)- 
valent and has D + 1 colorable edges. Here, instead, we 
choose to define a _D-colex A in terms of its dual lattice 
A. From now on D is an integer, _D > 0. 

Definition 1 A lattice A is a D-colex if its dual lattice 
A is a simplicial lattice with (D + 1)- colored vertices that 
forms a closed D-manifold. 

Recall that an n-simplex is the convex hull of a set of 
(n -|- 1) affinely independent points. In a simplicial D- 
dimensional lattice all the n-dimensional cells are home- 
omorphic to n-simplices. A closed manifold is a compact 
manifold without boundaries. For example, for _D = 1 it 



must be a circle, and for D = 2 [9| it could be a sphere, 
a projective plain, a torus, etc. That the vertices of A 
are {D + l)-colored means that we have attached labels 
to vertices from a set Q with D + 1 elements or colors, in 
such a way that two vertices connected by a link never 
share the same color. At first sight it may seem unclear 
whether a D-colex can be constructed in arbitrary man- 
ifolds, due to the colorability constraints. However, in 
[m a construction is given to obtain Z?-colexes from ar- 
bitrary lattices. 

We will denote the set of all n-simplices in A as A„, 
ri = 0, . . . , Z?. The coloring of vertices induces a coloring 
of 71-simplices. Namely, given a subset of colors q C Q 
with n + 1 elements, we define the subset A^ C A„ of 
those n-simplices such that their n+1 vertices are labeled 
with the elements of g. A coloring for the cells of the D- 
colex A is also induced. Let us write A* for the dual 
{D — ri)-simplex of a n-cell A and S* for the dual {D — n)- 
ccll of a ri-simplex S. Then if 5 is a q-ccll we say that 
(5* is a g cell, with q the complement of q in Q. As with 
simplices, the set of all n-cells of A is A„ and the set of 
q-cells is Ag C A„, q a Q, \q\ = n. Note in particular 
that vertices are colorless. Indeed, the colorability of 
vertices in A has its own content: vertices are bicolorable 
— that is, the graph of A is bipartite — if and only if the 
manifold is oricntable. In general A may be composed 
of m connected components A", a — 1, . . . ,m. We will 
regard each of these components as a '(Z?-|-l)-ceir, setting 
Aq -Ad+1 :={A"}. 

The previous notation for A and A can be extended to 
individual simplices and cells. Namely, (5„ C A„ will de- 
note the set of n-simplices that are part of an m-simplex 
6 E A,„, D > m > n > 0. For example, if 5 g Ai is a 
link then Sq = {e, 77} with e, rj the vertices in the ends 
of S. Similarly, A„ C A„ denotes the set of n-cells that 
are part of an m-cell A € A^, to > n, and Xq C A^ de- 
notes the set of q-cells that are part of a q'-cell A G Agi , 
qCq' CQ. 

An important feature of D-colexcs is that all their d- 
cells, < d < D, are {d— l)-colexes themselves. This can 
be easily checked in the dual simplicial lattice A. To this 
end, we first choose a vertex 5 G Aq. If we keep only those 
simplices that meet at S, and reduce their dimension in 
one by erasing the vertex 6, it is easy to check that the 
resulting simplicial lattice A' is dual to a (Z? — l)-colex 
A' with the topology of a (D — l)-sphere. This is the 
colex that correspond to the IP-cell 6* . The process can 
then be continued recursively to get the colexes attached 
to cells of smaller dimension. Thus, there is a one-to- 
one correspondence between subcolexes and cells. Each 
d-cell X £ Ad has as its boundary a (d — l)-subcolex. If 
A is q-colored, q = {(?;} C Q, then its {d — l)-colex is 
composed of gi-links, i ^ 1, . . . , d. 

Given a coloring r and a cell X E Aq, q C r C Q, 
there exists a unique r-cell, denoted A*" G A^, such that 
A € {X'^)q- It is natural to define the intersection of two 
cells fi G Aq, ly E Ar, q,r C Q, as the collection of disjoint 
{q n r)-cells v :— ^Jiqnr H Vqnr- When ^,C\v\s. empty 
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or has a single element for any cells ^ and v, we will say 
that A is nice. In a nice ZJ-colex, it fj, € Ag and X € ^ig 
with s d q d Q, there exists for any r <Z Q with s ^ qC]r 
a unique v G A,, with yu n = {A}, namely v ~ \^ . 



2. Color codes 

From a D-colex A we can obtain several color codes, 
labeled with an integer d such that Q < d < D. As we 
will see, the d-th code is dual to the J-th code. The 
codes with d = d are self-dual, something that turns out 
to have interesting computational consequences. 

To construct the c?-th color code from a ZJ-colex A, we 
first attach a qubit to each vertex v G Aq. As generators 
of the stabilizer S we choose certain {d+ l)-cell and {d + 
l)-cell operators, as we explain next. For any n-cell A € 
An, < n < D + 1, we introduce the X and Z cell 
operators 

:= n ■= n (^0) 

An important property of cell operators is that for any 
integers 1 < n^m < D + I with n + m > D + 1, and for 
any cells € A„, v £ A,„ we have [Tl| 

[X^,Z,]^0. (81) 

Also, for any g-cell A, q C Q, and for any coloring q' C q, 
we have the relations [ll| 

= n ^A.' = n (82) 

Definition 2 The d-th color code of A, < d < D, is 

the stabilizer code C with stabilizer S generated by the 
operators of the form X\, A g A^+i, and Zy , A' G A^^^. 

From this definition it is clear that dual codes are related 
through the exchange of X and Z operators, and thus are 
equal up to a change of basis. 

From the relations (|82p it follows [ll| that the number 
of encoded qubits is 




where = kg is the d-th Bctti number of the manifold. 
The origin of this particular value will be clarified in the 
following sections. 

Note that color codes with d = or d = D are not 
interesting quantum codes, as they are repetition 'classi- 
cal' codes, which only correct either X or Z errors and 
encode a qubit per connected component of A. Because 
of this, below we will mainly focus on color codes with 
0<d<D. 



3. A simplicial homology 

In this section we will work only in the dual lattice A, 
since it is more suited to explain the underlying homo- 
logical structure of color codes. We attach to each set of 
ri-simplices A„, Q < n < D, the abelian group C„ of Z2 
chains, together with a trivial group C_i = {0}, and to 
each n = 0, . . . , D a boundary operator 

Cd ^ C^-i (84) 

defined for n > by 

where (5 e Ajj. The exact sequence from which the c?-th 
color code is obtained is 

Cd-1 — * — ^ Cd-1, (86) 

where < d < D. It gives rise to the homology group 

Hd := Za/Bd. (87) 

It follows from ([55]) that the homology groups H^, Hd 
are isomorphic to the the sum of {^)hd copies of Z2. In 
other words, we have 




The origin of the binomial coeffiecient is related to color 
and will be explained in the next section. 

4- Brane-nets 

We will now give a geometric interpretation to the D- 
simplex cycles discussed in the previous section in terms 
of branes. For any n-cell A G A„ we define cx & Cd as 

CA := E (89) 

veXo 

For < n < D, we will now introduce boundary opera- 
tors that take n-cells to (n — l)-chains in a 'natural' way. 
For A e A„ we set 

dX:=dnCx. (90) 

We say that these definitions are natural because dX, 
as a set of simplices glued together, is homeomorphic 
to the boundary of A in the usual sense. Moreover, if 
A is a g-cell, q <Z Q, then dX only contains g-simplices. 
Thus, we can regard n-cells as pieces of n-branes. Indeed, 
coloring is important here, and we should regard a g-cell, 
q C Q, as a. piece of a g-brane. As we will see next, such 
pieces can be put together to form closed q-branes, that 
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is, branes without boundaries in the sense of the d just 
defined, which is just a reinterpretation of i9„. 

Let us fix a particular color subset q with |g| = n and 
introduce the g-colorcd chain subgroup C|) , the subgroup 
of Cd generated by all c\ chains with A G Ag. That is, 
the chains in C|) represent the g-branes just discussed. 
We may then set Zq := Z„ n C^, Bq B„ n C^, and 
Hq := Zq/Bq with the net result 

Hq ^ Hn. (91) 

Thus, the homology of g-branes for a fixed q is just the 
usual brane homology. This result is shown in [llj by 
constructing from the colex A a g-reduced lattice in which 
all n-cells correspond to g-cells of A. 

In summary, the cycles in Zq can be regarded as closed 
g-branes, which arc also cycles in Z„. In general, cycles 
in Zn will combine several pieces of such g-branes of dif- 
ferent color, forming n-brane-nets [Tl] |. In this sense, iJ„ 
is a n-branc-net homology group, as opposed to the usual 
Hn homology group of n-branes. This discussion already 
suggests that iJ„ contains several copies of iJ„, but how 
many of them? To answer this, let us consider color sub- 
sets qi C Q, i = 1, ... I? + 2 — n, with l^j = n and 
\qi r\ qj\ = n — 1 ior i ^ j . It turns out [ll|] that if Zi, 
i = 1, . . . D + l~n, are g^-branes, all of them homologous 
in the usual sense of Hn homology, then 

^2,-0 (92) 

i 

up to Hn homology, which shows that only (^) color com- 
binations are independent, which finally explains ()87p . 

In a d-th color code the relevant homology groups are 
-ffj and Hd, with the corresponding d- and o?-brane-nets, 
which are mapped to operators using (|64p . The com- 
mutation rules of closed brane-net operators are purely 
topological. Namely, given a closed g-brane-net z ^ Zq, 
\q\ = d, and a closed g'-brane-net z' e Zqt, |g'| d, the 
operators and Zz' anticommute if z and 2' intersect 
and odd number of times and g H g' = 0. Otherwise they 
commute. 



5. Branyon-nets 

For d = 1,...,D, the Hamiltonian model ([75|) of a 
d-th color code contains as in toric codes two kinds 
of brane- like excitations, (d — l)-branyons and {d — 1)- 
branyons. These are created with the brane-net oper- 
ators discussed in the previous sections, in particular 
with open branc-nets since the excitations appear at their 
boundary. Then, unlike in toric codes, branyons can form 
nets and we must talk about branyon-nets. As branes in 
the previous sections, branyons can be labeled with col- 
ors. Namely, a g-branyon is the boundary of a g-branc. 

We will only discuss (d— l)-branyons, which are related 
to violations of the stabilizer conditions at {d + l)-cclls. 



The case of {d — l)-branyons is analogous. If A S ^d+i 
is a g-cell, we say that the corresponding excitation is 
a (Q — g)-branyon excitation. We know that branyon 
excitation must be arranged in closed structures. This 
fact is reflected [ll| in a 'conservation law' at {d + 2) 
cells. Namely, for A a g-cell, |g| — d + 2, and for any 
r, s G q, \r\ = |s| = d + 1 we have 

n = n (93) 

which is just a special case of ([5^ . This implies, for 
example, that at a given (d -I- 2) cell there cannot exist a 
single excited {d+ l)-cell. 

In the previous section we have seen that only (^) 
color combinations are enough to label an independent 
set of d-brane operators. We will now show that we can 
locally transform a set of branyon excitations into a new 
set that only contains branyons from (^) independent 
colors. By 'locally' wc mean that it is enough to act 
in the neighborhood of the undesired excitations. We 
start choosing a color go € Q- We want to show that 
we can eliminate all g-branyon excitations with go G g 
or, to put it in a different way, all excitations located at 
cells A € ^d+i of color g with go ^ g. The total set of 
excitations corresponds to the boundary of some c £ Co, 
that is, they can be switched off with Xc- Notice that all 
the g-cells with go ^ g that touch a (Q — go)-subcolex are 
part of it. Then the restriction of Xc to those {Q — go)- 
subcolexes that contain the undesired excitations will do 
the job. 

Once we have reduce the branyonic excitations to the 
(■|^) color combinations, each color may be treated in- 
dependently as in the previous section. That is, the 
branyons of each color form boundary (d — l)-branes. 

Finally, we note that, as in the case of generalized toric 
codes, one can derive |ll| from the commutation rules for 
closed brane operators the topological ±1 phases appear- 
ing in branyon processes. 



6. Self- correcting models 

Regarding self-correction, color code models behave 
very much like toric code models. That is, for d > 1, 
the models with Hamiltonian (|75|) for the d-th code pro- 
tect Z-type encoded operators, and for d < Z? — 1 they 
protect X-type operators. Instead of repeating the anal- 
ysis already done for toric codes, we will just mention the 
new features that appear in color codes. In particular, 
wc consider the correction of A"-type errors. 

The first difference is the lattice, but any periodic D- 
colex will do, and the existence of periodic lattices follows 
from the procedure given in [llj to construct D-colexes 
from arbitrary lattices. The second difference is that now 
we deal with branyon-nets instead of branyons, but this 
does not make such a big difference due to the discussion 
in the previous sections. In this regard, we first need a 
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notion of connectedness of {d + l)-cells. This is already 
implicit, for example, in (|93p : we construct a graph with 
l)-cells as its nodes in which nodes are connected by 
a link when the corresponding cells are part of the same 
{d + 2)-cell (notice that d > 2). As for the corr function, 
we transform each connected set of excitations using the 
procedure of the previous section into a collection of {d — 
1) boundaries of (^) colors, so that each color can be 
independently treated exactly as in a toric code. 

F. The transversality properties of color codes 

Color codes can have very special transversality prop- 
erties when the lattice and its topology are suitably 
tuned. In this section we will treat this issue in detail. 
To this end, we start explaining a simple way to consider 
D-colexes in manifolds with boundary. We then describe 
a particular topology for the manifold: in the shape of 
a D-simplex in which each face has different properties. 
Finally, we analyze the transversality properties of the 
resulting simplicial codes. 

1. Boundaries 

Up to now we have only considered colexes A in closed 
D-manifolds. However, richer color codes can be ob- 
tained by introducing suitable boundaries in the manifold 

[l3l . The easiest way to describe such boundaries is as 
collections of missing cells in the lattice. For a missing 
cell in A we really mean a missing generator of the stabi- 
lizer. When the removed generator is independent from 
the remaining ones, the number of encoded qubits k gets 
increased (jH). This can be understood also in homologi- 
cal terms. Consider a g-brane c G Co with a boundary, 
\q\ ~ d. Given an encoded state the state Xc\ip) will 
have stabilizer violations at those {d+ l)-cells that form 
the boundary of c. But if we eliminate the generators 
that correspond to these {d+ l)-cells then Xc will belong 
to the stabilizer and thus we should consider c a closed 
g-brane. This way, as cells are eliminated we are intro- 
ducing in the manifold areas in which suitable branes can 
have a boundary and still be closed. It is as if we had re- 
moved a piece of manifold in these aereas, thus creating 
a boundary. Such an introduction of a boundary does 
not change homology only in the way just discussed. In 
addition, when we remove a g-cell the corresponding gen- 
erator, which can be regarded as a d-brane, transforms 
into an encoded operator in Af — S. Thus, this d-brane 
should no longer be considered a boundary, but rather a 
nontrivial cycle. 

Although for our purposes it will be much more con- 
venient to describe boundaries as big missing cells in the 
lattice, this would not be convenient at all in an imple- 
mentation of the model. Fortunately it is pos sible to 
adopt a different approach, as discussed in [36|[in the 
context of 2D toric code models and also in [33 in a 



more general framework. In this alternative approach 
boundaries appear at interfaces between the topologi- 
cally ordered phase and other phases in which some of 
the branyons are condensated. 

2. Simplicial color codes 

There exists a very special class of color codes which 
live in simplex-shaped lattices. For D = 2 we have tri- 
angular codes 0, which bear this name because their 
lattice is a triangle: each side of the triangle corresponds 
to a different kind of boundary. For I? = 3 we have 
tetrahedral codes [To| . and more generally for arbitrary 
D we may consider simplicial codes. Instead of giving 
a detailed description of the properties of each of the 
boundaries in a simplicial code, it is just easier to give 
a constructive procedure to obtain such codes, as we do 
next. 

Suppose that we have at our disposal a Z?-colex A with 
the topology of a D-sphere, D > 2. Since spheres are 
homologically trivial for D > 2, A yields trivial color 
codes with fc = 0, no encoded qubits. Let us now choose 
a vertex of A and proceed to remove the corresponding 
physical qubit and all the {d + 1)- and (d -I- l)-cells that 
meet at this vertex. The resulting lattice looks like a D- 
simplex, because it is a D-haXl with its boundary {D — 1)- 
sphere divided in the D + 1 regions that where in contact 
with each of the D + 1 removed Z?-cells. It can be checked 
using (|82p that only two independent generators of S are 
eliminated in the process, one of X type and the other of 
Z type. Thus all color codes derived from the new 
lattice A' encode exactly one qubit. Since the original 
lattice must had an even number of qubits ^ic|| , it follows 
that the operators 

X:=A®", Z:=Z®", (94) 

with n the total number of remaining qubits, are the en- 
coded X and Z operators: they belong to the normalizer 
and 

{A,Z} = 0, (95) 

because the original number of qubits is even (|82p . 

It is easy to construct the smallest example of a simpli- 
cial color code for each dimension D > 2. We start with 
a (I?-|- l)-hypercube, in which we color parallel links with 
the same color to get [ll| a D-co\ex with the topology of 
a D-sphere. Then as explained above we eliminate one 
of its vertices and all neighboring cells. As a result we 
get a color code with 2^+^ — 1 physical qubits in which 
all n-cells have 2" vertices. 



3. Transversal gates 

We will now consider the tranversality properties of 
color codes. As it turns out, simplicial colexes give the 
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most suitable color codes for the transversal implementa- 
tion of quantum gates. Thus, we restrict our discussion 
to them. 

First, it is clear that the transversal implementation 
of the logical X and Z operators is always possible, 
since they take the form (|94|) . The same reason makes 
transversal measurements of X and Z are always possi- 
ble. Similarly, the two-qubit controlled not (CNot) gate 



C 



X 



\1){1\®X 



(96) 



can always be implemented transversally on pairs of equal 
simplicial codes by applying it to each pair of equivalent 
qubits. The one-qubit Hadamard gate 



V2 



1 1 

1 -1 



(97) 



is more restrictive. As it exchanges X and Z operators, 
it can only be apphed to dual codes, those with d = d. 
Again, it can be obtained by applying a Hadamard to 
each of the physical qubits independently. 

We want to consider in general the family of one-qubit 
phase gates 



Rk 



|0)(0| 



-exp(J)|l)(l| 



(98) 



where fc is a non-negative integer. The case fc = is 
just the X gate. The case A: 1 is important since 
the gates {Cx , H, Ri} generate the Clifford group. Still 
more interesting is the case k = 2, because Cx and i?2, 
together with X and Z intializations and measurements, 
arc enough for universal quantum computation[38j . To 
study in detail when these gates can be implemented 
transversally on simplicial color codes, we will introduce 
next the notion of (j, A:)-good colexes. 

From now on, j, k are integers with < j < D, k > 0. 



Definition 3 A D-colex A is {j,k)-good if 



Vc e B 



J' 



= mod 2 



fe+i 



(99) 



Note that all colexes are {j, 0)-good. From the first ex- 
pression for encoded states in (|7ip and the fact that sim- 
plicial colexes have an odd number of physical qubits it 
follows that simplicial color codes obtained from (d, fc)- 
good colexes allow the transversal implementation of Rk 
gates. More exactly, if we apply Rk to each physical qubit 
independently we will get a logical Rf. gate, with an odd 
power s. But since s and 2*^+^ are coprime, there exists 



an integer r such that R^^ — R^. 



Rk- 



The problem with definition [3] is that it is not really 
operative. But we can give a simpler characterization of 
(j, fc)-goodness. For its proof, see appendix IB] 

Theorem 4 A nice D-colex A is (j, k)-good for some 
k > if and only if j' := 2j — D > 0, all j-subcolexes are 
(j', k — l)-good and 



VA e A 



|Ao| = mod 2*^+1, 



(100) 



Moreover, the <= implication holds even if A is not nice. 

Since (0, fc)-good colexes only exist for /c = 0, by re- 
peatedly applying theorem [5] it follows that a (j, fc)-good 
and nice _D-colex may exist only if 



j> 



k + 1 



-D. 



(101) 



Several conclusions may be derived from this result. For 
example, that dual codes only allow the transversal im- 
plementation of the i?i gate. It could not be other way 
else, since if a dual code allowed the transversal imple- 
mentation of i?2 wc would get a stabilizer code with a 
universal set of transversal unitary gates, something im- 
possible [1^. For fixed J, there is a minimal dimension 
D for which (d, fc)-good nice D-colexes exist, namely 



D > {k + l)d 



(102) 



In particular, setting d = 1 it follows that the minimal 
dimension D for which Rk may be transversally imple- 
mcntable in a color code is D = fc + 1 . 

Once the bound (|102p is established, it only rests to 
show that it is saturated by some lattice. But the sim- 
plest examples of simplicial lattices given in the previous 
section do indeed saturate the bound. The next question 
is whether periodic lattices that saturate it exist. This is 
indeed the case for D ~ 2,3, and we will assume that it 
is also true for higher dimensions. 



VIII. SELF-CORRECTING QUANTUM 
COMPUTERS 



A. Initialization from thermal states 

In the previous sections we have studied how quantum 
memories resilient to thermal noise can be obtained from 
local quantum models in 4D and above, both from toric 
codes and color codes. However, a quantum memory is 
of no much use if we cannot initialize it. We will now 
show that it is possible to initialize the self-correcting 
topological quantum memories as long as we can control 
the local quantum Hamiltonian. In particular, we will 
consider toric and color code models with d > 1, so that 
the dressed version of Z-type encoded operators TV are 
protected, and give a procedure to initialize all such oper- 
ators to iV = 1 with a probability of failure exponentially 
small in the system size. For D = d = 2, this is just the 
usual initialization of a 2D ferromagnet with a vertical 
magnetic field. 

Our starting point is a modification of the Hamiltonian 
((75)) . We introduce single qubit Zeeman terms and some 
couplings to reflect the fact that now the strength of the 
different terms can be controlled: 

I m 

H = -9x Xp^ ~9zY. Z,^ - go ^ Z^, (103) 
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where s runs over all qubits and gx,gz,go > 0. Be- 
fore considering the details, let us sketch the initializa- 
tion procedure. We start from a thermalized state with 
gz = go = ^ E^nd gx = 0. Then we switch off the Zeeman 
terms in some time Ti. Finally, we 'carefully' switch on 
the X-type stabilizers in some time T2 till 51 = 1. No- 
tice that the Hamiltonian is exactly solvable at all times, 
since all the terms with non-zero coupling at a given time 
commute. Moreover, all the terms commute with Z-type 
encoded operators, precisely the ones that we want to 
initialize with eigenvalue 1. 

We first want to show that, for temperatures 1/(3 below 
a critical value, the initial thermalized state satisfies, for 
any encoded operator iV e A/" n Vz and some < > 0, 

1- (iVdr),3 < e-*^. (104) 

Recall that gx = Q and .gz = 5o = 1- Both for toric and 
color codes we already have a graph of the excitations 
that correspond to gz terms, where nodes are labeled 
by the qi elements of the lattice. We enlarge this graph 
by adding a node per physical qubit, to represent the 
excitations that correspond to the Zeeman g^ terms. We 
also add a link from the node of the physical qubit s to 
the node of qi if {X^, Z^^} = 0, and a link between the 
nodes of two physical qubits s and s' if they are both 
linked to some qi. For example, in a toric self-correcting 
code with d ~ 4 the Qi nodes correspond to edges and the 
s nodes to plaqucttes, so that the links of the graph just 
represent the usual notion of connectedness for plaquettes 
and edges in a lattice. The Zeeman term is responsible for 
the tension of membranes built of plaquettes, while the 
Kitaev term for the tension of the loops, i.e., boundaries 
of the membranes. 

It is then clear that any connected set of excitations 
can be independently switched off with a suitable opera- 
tor in Vx- We can then repeat the reasoning of section 
IVI CI to show that those excitation configurations which 
contain connected sets of size bigger than, say, L/2 are 
exponentially suppressed with the size of L. Therefore, it 
suffices to observe that the rest of configurations belong 
to the sector with iVdr = 1 for any iV G A/" n Vz, as we 
detail next using the language of section lVII CI First, the 
eigenstates of the Hamiltonian may be written as \a)o for 
some a £ Ai. The corresponding set of excitations cor- 
responds both to the chain a and its boundary dia £ Aq. 
But corr(S'2 n da) = Xa' for some a' G Ai such that both 
a and a' lie in the same hypercube of length L/2, and 
thus XaXa' e S. The result follows. 

Once we have shown that the initial thermal state is 
in the desired sector, it is enough to check that this con- 
tinues to be true as we switch off and on the interactions 
as indicated. This follows from the fact that gz = 1 at 
all times, as long as we keep the temperature of the sys- 
tem stable. In this regard, the process of switching on 
gx could be a problem since at start we have arbitrarily 
many violations of the Xp. stabilizers and if gx raises 
too fast we will put a lot of energy in the system. To 
avoid this, we require that the switching on is done slow 



enough for the system to progressively thermalize. We 
assume that the thermalization time is size independent. 

B. Operations 

From the discussion in section lVII F Sl it follows that 6D 
color codes are enough to have a set of universal transver- 
sal operations. Local CNots require an additional dimen- 
sions, so that we need 7 spatial dimensions. 

1. Measurements 

The fact that we are dealing with CSS-like codes im- 
plies that we may perform transversal X and Z mea- 
surements. Typically, measuring a physical qubit in the 
Z basis requires an interaction Z (g) Happ with a mesur- 
ing apparatus. In this regard, during the measurement 
we can switch off the X-type stabilizers, which do not 
commute with the measuring Hamiltonian. Switching off 
this part of the protecting Hamiltonian is not a source of 
problems since the corresponding errors do not affect the 
measurement. 



2. Transversal gates 

In order to perform transversal gates we have to tem- 
porarily modify the Hamiltonian of the system. This 
means adding new terms that will not commute with the 
protecting ones and thus make the analysis difficult. In 
any event, when a gate is half-way of being completed the 
system will in general be in a state with a high energy 
with respect to the protecting Hamiltonian. This could 
easily be a source of problems. 

To circumvent them, we propose to switch off the pro- 
tective Hamiltonian while performing the gates. We will 
argue that the gate can be applied with high fidelity to 
the logical qubits, as long as it can be done fast enough 
(of course at some constant speed not depending on the 
required fidelity of the gate). First, notice that we can 
represent the faulty gate as an ideal gate preceded by 
some noise operation, so that it is enough to consider 
a situation with a trivial gate (identity), i.e., qubits are 
subject solely to the interaction with the environment. 
Thus, the process can be seen as a noisy channel, so that 
we can use the results of section IVI Dl which tell us that 
the quantum memory can resist a small amount of noise. 
Indeed, the results in that section where obtained for a 
depolarizing channel acting independently on each qubit, 
but they can be generalized to more complicated noises 
as long as correlations are local. 

However, this is not really the end of the story. Even 
if we could perform transversal gates perfectly and in an 
instant, the gates themselves will increase the number of 
excitations. In the case of single-qubit gates we know 
that excitations cannot spread much, so this is not a 
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big problem. The case of CNot gates is worse, because 
they copy errors (and thus excitations) from one code 
to another. It is intuitively clear, however, that as long 
as excitations are scarce enough this will not affect the 
encoded qubits. A rigorous approach to this problem 
requires using similar arguments to those in sections I VI CI 
and IVI Dl but considering two copies of the Tl lattices 
together. 

In summary, the application of a transversal gate will 
put some amount of energy per volume unit in the sys- 
tem. As long as the temperature is low enough and the 
energy not too much, this will not produce an error. Of 
course, we have to wait some time between gates, so that 
the injected energy can be dissipated. 



IX. CONCLUSIONS 

We have constructed a model for a self-correcting 
quantum computer based on topological stabilizer codes. 
The model, to be local, requires 7 spatial dimensions. 
Our model consists of 1) a lattice of qubits exhibiting 
proper topological and geometrical properties so that 2) 
a code on this lattice allows a universal set of transver- 
sal gates (including measurements) and 3) a protect- 
ing Hamiltonian, with controllable coupling constants. 
We have described the procedures of initialization and 
measurement of the quantum computer, and how gates 
should be performed. For the latter, we propose to tem- 
porarily switch off the protecting Hamiltonian so that we 
can avoid problems due to the combined effect of the pro- 
tecting Hamiltonian and the one needed to perform the 
gates. 

We have then obtained two rigorous results concerning 
the self-protection of such a quantum computer. First, 
we have shown within the weak coupling limit regime 
that the memory is stable, generalizing thereby results of 
. Second, we have shown that if the memory is exposed 
to a local noise for a short time the quantum information 
is preserved. Our further results are not rigorous but are 
physically well motivated. These comprise the initializa- 
tion scheme and a reasoning showing that we can perform 
gates on protected qubits with high fidelity. 

We would like to emphasize that as the formulation 
of a self-correcting quantum computer is a challenge, we 
decided to relax the number of spatial dimensions D in 
which our proposals are formulated. The point is to ad- 
dress this problem as a fundamental issue regardless of 
the dimensionality. Once we have seen that this is pos- 
sible under the circumstances considered in this paper, 
then we may reflect on possible ways to lower the dimen- 
sionality D where the lattices are embedded. To this end, 
we should consider our models as simplified models as far 
as the types of degrees of freedom is concerned. We have 
only used spin 1/2 or qubits. It is an open question to 
see whether by including additional degrees of freedom 
like fcrmions or bosons, we could induce interactions that 
may lower the dimension in which self-correcting quan- 



tum computers are possible • A similar possibility is to 
consider models with generalized spins and discrete non- 
abelian gauge symmetry groups [40] which are found to 
posses string tension and confining properties. Generally 
speaking, there are instances of quantum theories whose 
critical dimensions for being stable can be lowered by 
introducing additional types of degrees of freedom or a 
combination thereof. This could also be the case here and 
this possibility we leave it open, but with the new insight 
gained here, we see that at least with simpler models it 
is possible to achieve the property of self-correctedness. 

Finally, we should discuss the apparent contradiction 
between our results and and those in ref. [l3|. In that 
paper Robert Alicki argues that a quantum computer 
could be used to draw work from single heat bath. To 
this end he assumes that a protected bit cannot be used 
to draw work from a heat bath, i.e. it docs not follow 
the Landauer principle. The author implicitly does not 
allow to switch off the protective Hamiltonian. Indeed, 
when the protective Hamiltonian is all the time on, the 
protected qubit does not interact with the reservoir, and 
there is no way to use it as an entropy sink in order to 
draw work from a single heat bath [4l|. On the other 
hand, if the protective Hamiltonian can be switched off, 
the qubit can interact with the bath, and it can be used 
to draw work (see e.g. [i^). Now, in our scheme, we 
need to switch off part of the Hamiltonian during state 
initialization, and therefore the conclusions of [l3| do not 
apply. It is interesting that the need of switching off and 
on the protective Hamiltonian for the quantum computer 
have been obtained from completely different points of 
view. Alicki obtains it from an elegant thermodynamical 
considerations, while we need the switchablc protective 
Hamiltonian in order to avoid non-commuting terms both 
at initialization and at gate application. Thus, in our ap- 
proach, the need for switchable protection was a matter 
of convenience (otherwise, we simply would not be able 
to present a plausible scheme), while Alicki's argument 
support the thesis that this property is actually necessary 
in such a quantum computer. 
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APPENDIX A: A BOUND FOR CONNECTED 
COLLECTIONS 

We want to show that, in a graph with at most fi 
links meeting at each node, the number of connected col- 
lections of nodes of size k that contain a give node is 
bounded by e^'' for a suitable r that only depends on ii. 
We first notice that it is enough to consider a /i-valent 
tree (infinite), because this is the worst case. Then a 
connected collection of nodes is a subtree. The idea of 
the proof is to express a subtrees T as a word in a 
suitable alphabet with 2^* letters. For a given node, we 
will suppose that its neighboring nodes are ordered. 

First, we know that certain node n is part of the sub- 
tree. There are fi nodes connected to n. The possible 
subsets of these nodes are 2^. The first letter of W tells 
us which is the subset contained in T. Then we jump to 
the first node n' of these collection, and add its neigh- 
boring nodes as dictated by the second letter. Then we 
jmirp to the first of these new nodes, and so on. When 
there are no more nodes to add at a given node, we jump 
to the last node that was added. When there are no 
more nodes left, the construction ends. Since a tree with 
k nodes is expressed with a word of k letters, the final 
result is r = //log 2. 

APPENDIX B: PROOF OF THEOREM H 

First we need some extra notions for those cases in 
which j' > 0. Recall that each A E ^j+i is a j-colex A'. 



We introduce a group homomorphism hx : Cd — > Cj 
that restricts c € Cjj to Cj , the j-simplex chain group 
of A. In particular for c = X^ugv ^* "^ith C Aq we 
set hx{c) := Y^vev^* '^i*^ V ^ V Aq. This homo- 
morphism maps Bj to _Bj, . To check this, it suffices to 
show that for any /i e Aj+i we have /ia(c^) = Siyew 
for some N C A^'+i. But if /.( n A = {v^} we just set 

When A is nice, we can in addition define a group 
homomorphism : B'^, — > Bj such that for every c' e 

B'j, we have hx{h^{c')) = c'. For each v e A^v we set 
h'^{c'^) := c,jq, where q ~ Q ~ [r — s) with s,r C Q such 
that 1/ e As, A e A^. Here c'^ is meant to denote hx{cu). 

We can now prove each of the implication directions: 

/ First, for any A G A_,_|_i we have cx € Bj, so that 
|Ao| = |ca| = mod 2''+^. Suppose that f < 0. Then 
for any two colorings q,r d Q with \q\ = \r\ — j + 1, 
|(7 n r| =0 and for any u G Ag we get, using the niceness 
of A, that |c„5 + c^H = Icu-J \ + \cy^ \ - 2 = -2 mod 2'=+^, 
a contradiction. Finally, suppose that A G A^j+i is not 
(j',fc — l)-good as a j-subcolex. Then, with the same 
notation as above, there exists a chain c' G Bj, with 
|c'| ^ mod 2*^. For c = h^(c') we have by assumption 
|c| = mod 2'=+!, but then \c + cx\ = |c| + |ca| - 2|c'| ^ 
mod 2'""+^. Since c, ca G Bj, this is a contradiction. 

<;= / It suffices to show that for any c G Bj with |c| = 2 
mod 2*"'+^ and for any A G A^+i we have |c + ca| = 
mod 2*-'. But for c' = hx{c) we get |c -I- cx\ — |c| + |ca| — 
2\c'\ and since |c'| = mod 2^ the result follows. 
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